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^ . Abstract 

in . 

^^ I One of the necessary steps in constructing high precision option of KKMC was to 

\^ • instah the double bremsstrahlung matrix element for the process e~^e~ — > fgj^e into 

^^ , the scheme of Coherent Exclusive Exponentiation. The process is also interesting 



o 



Ph: 



because of gauge cancellation of contributions for photon emission from incoming 
fermion lines and t-channel W. The QED U(l) gauge properties require terms of the 
triple and quatric gauge couplings to be taken into considerations as well. Thanks 
, to expansion starting from the approximation of contact interaction, good example 

^' to study the internal structure of the amplitude is available. 

r^ • In the developed scheme, natural separation of the complete amplitude into 

gauge invariant parts is straightforward. Each part has well defined physical inter- 
is^ I pretation, which after partial integration over phase space provides terms: infrared 

j_j ' singular, leading log, next-to-leading-log, etc. Contributions related to triple and 

^ quatric gauge coupling of W (extracted with the help of expansion around contact 

VF-interaction) , have been ordered as well. The separation is also helpful, to define 
extrapolation/reduction procedure of CEEX exponentiation for the i^e channel. 
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1 Introduction 

Higher order radiative corrections are usually necessary to obtain from Standard Model 
high precision results for phenomenologically important quantities. The techniques of 
direct calculations, lead to expressions, of often hundreds, thousands or even million's of 
terms. These are difficult to control analytically and/or numerically. This is worrisome, 
because to obtain phenomenologically sound results third order effects are mandatory, 
see e.g. [1]. This is clearly outside the reach of presently available methods of direct 
perturbative calculations. There is no doubt, that resummation is necessary of at least 
some contributions from orders higher than the second one. 

In case of electroweak processes at LEP techniques based on exclusive exponentiation 
of QED effects turned out to be powerful and enabled high precision predictions for a 
wide game of processes, such as Bhabha scattering, production of heavy bosons W or 
Z and lepton pairs. The underlying method, originating from pioneering work of Yennie 
Frautschi and Suura [2] turned out to be realizable [3-7] in practice, thanks to acumulated 
experience and ever increasing computer power. 

One of the necessary elements in approach based on exponentiation is rigorous study 
of matrix elements obtained from perturbative calculation. In fact it is not enough to 
calculate predictions at as high order of perturbation expansion as possible, but also to 
carefully separate results into infrared singular and remaining finite parts. Thanks to 
the properties of QED, results of explicit perturbative calculations are not necessary to 
obtain singular and leading terms of every order. These leading parts of amplitudes can 
be combined with the phase space into the module of the low level Monte Carlo generator 
or in general into multi-dimensional distribution which can be understood as lowest order 
of improved perturbative expansion. Later, finite parts of the matrix elements can be 
added order by order. In case of Monte Carlo algorithms it can be done with the help 
of the correcting weight, which can be shown to be positive and bound from above. 
Details of such a rigorous scheme can be found in refs. [7,8]. It improves significantly 
the convergence of perturbative expansion; final states with arbitrary number of photons 
are present already at the lowest level of expansion. This allows for predictions with 
realistic experimental cut-offs included. The solution based on separation performed at 
the amplitude level, is specially useful. It opens the way for easy implementation of all 
sort of interferences, also convergence of expansion is particularly fast in this case. The 
underlying exponentiation scheme is called Coherent Exclusive Exponentiation (CEEX). 

The following point is of practical importance. In case of exponentiation, already 
at the lowest order of expansion, configurations with multiple real photons are present. 
That is why, it may happen, that for particular event, there is more explicit photons in 
final state, than in expression directly available from standard perturbative expansion. 
Reduction/extrapolation methods are then necessary. We will not elaborate much on 
theoretical aspects of this point here. However let us stress that if sufficiently high order 
of perturbation expansion is available, dependence on the choice of reduction procedure 
or extrapolation is dropping out. Particularly bad choice may nonetheless degrade con- 
vergence of expansion. It is thus of the importance to provide results of perturbative 



calculations in a form as convenient for extrapolation procedure as possible. Comparisons 
of amplitudes calculated at different orders of perturbative expansion can provide a useful 
hint. 

In the present paper, we heavily rely on ref. [9]. We will assume certain level of 
familiarity of the reader with that reference. Also, let us note that spin amplitudes for 
the process e^e~ — > z/ei^e77 are well defined within the Standard Model and known since 
long, see e.g. [10]. We could profit in our work from the ready to use computer codes such 
as [11] available for numerical cross checks of our results. 

Our paper is organized as follows. Section 2 is devoted to the case of single bremsstra- 
hlung e+e~ — *> t'e^'eT, some basic elements of spin amplitude techniques, useful in the 
more complex case of double bremsstrahlung are presented there as well. In section 3, we 
provide the main results, in particular we explicitly identify gauge invariant parts of the 
amplitudes. We stress points which will be useful for extrapolation schemes used in CEEX 
exponentiation as well. We keep our discussion, having in mind future applications in spin- 
amplitude automated programs. In section 4 we discuss issues related to extrapolation 
procedure in more detail. Finally, Section 5 closes the paper. 

2 Amplitude for one real photon and notation 

Let us start with the well-known and straightforward to calculate by any method 0{a) 
spin amplitude for the e+e^ -^ i^ei^el single-photon bremsstrahlung process, see fig d 
We will recall it nonetheless here, to define framework for our discussion. We will use 
conventions of refs. [8,12]. Let us recall here only the most important notations. The four 
momenta Pa, Pb, Pc, Pd, ^i denote respectively momenta of incoming electron, positron, 
outcoming neutrino, antineutrino and finally photon. The indices for the spin states for 
the fermions are denoted respectively as A^, Af,, Ac, A^ and for photon cti. The photon 
polarization vector is denoted as eo-i- The gauge transformation in our case reduces to 
the replacement eo-^ — >■ eo-i + a; ki (with arbitrary coefficient x) , there will be no external 
line bosons, and incoming fermions lead to the trivial phases only. With these notations 
the first-order matrix element^ obtained from the Feynman diagrams depicted in fig. ^ 



^Mi{/} (actO The subscripts 1 and {/} denote respectively, that the amplitudes are of the first order 
and are included as part of the initial state bremsstrahlung. This spurious notation is however convenient 
for the reader interested in ref. [9]. 
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Figure 1: The Feynman diagrams for e+e 
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can be written in a rather straightforward way: 
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or, equivalently: 
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where, the part of the amphtude, consisting of bosonic couphngs (g^ denote couphng 
constant of Z with fermion / and handedness A, in electric charge units), final state 
fermion spinors and boson propagators reads as 
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The final-state spinors are explicitly included, and Fierz transformation is applied for 
the part of W exchange. The W coupling constant reads 

y/2 sm &w 
Only for the W contribution, the superscripts xy in M{/} have the meaning, they define 
the momentum transfer in the W propagator n|^(X): for xy = ac the transfer ^ is 

^Transfers can be expressed also as ta = (pb — ki — pd)^ and % = {pa — fci — Pc)'^, this make difference 
if extrapolation procedures are used for the configurations off mass shell where Pa + Pb 7^ Pc + Pd + ki, 
otherwise Al^ = A4^ of course. 



ta = {Pa — PcY, for bd it is tf, = {pf, — pdY- If both are explicitly marked, then the 
expression 

M5f;^ = ze2G,^^(Gr''^)M] (6) 

is used. For that parts of formula (j21) W propagators are explicitly given. The notations 
TZz and IZw will be used later, see respectively formula (j^ and (jTfjl . 

Let us start now to rewrite expression (J2I). It is straightforward to notice that the 
first term M.^ can be split into soft IR parts proportional to (/^ ± m) and non-IR parts 
proportional to ^i. The non-IR parts are individually gauge invariant by construction. 
The soft part of A^°, with Z couplings only, is gauge invariant as well. 

Employing the completeness relations of eq. (A14) from ref. [8] we obtain the different 
form of 0: 

-''"-i{^} KxaJ - 2k p ^ ^^''^"^ ' ' U^iaJ ^ 2k Pi, ^ LAi,aipJ ^ LpbAj M] 

pa Pb 

eQe^ \^mrPbfcil, rr pifciPal ^^e V^ t^ Tpb fci fci] qi plPa] ^ ^ (7) 

+^i::;r 2^ ^ La,p J mf U -iaJ - t^t— i^ v U,-ip J ^ U aJ m 
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The terms Ai^ to Ai^ correspond to the last three lines'^ of eq. (^. These contributions 
are also IR-finite. In the next step let us remove the sum in the first two terms thanks to 
the diagonality of U and V ( ref. [8]). The matrices 03 are also defined in this reference. 
We obtain 

+ 4^ + rf;^ + {M' + M') 

„B' (Pkl\ _ ^Qe V^m \PbPa-\ ^ TJ \PaklVa^ 



2kiPa 
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„S" (Pki\ _ , (^Qe \^Y \PbkiPb] m \PbPal , „ 
^{1} [Xai) - + 2k^pf, ^ LAtaipJ -^ [pbXai M 

A' (pki\ _ , ^Qe V^a> rPbfci] rr ffcifciPa] 

2kxPa 2kipb 



(8) 



^The term M^ + M.^ + M^ originates from the WW^ vertex 

-ie [gt^vip - q)p + 9up{q - r)^ + g^p{r - p)^] 

where all momenta are outcoming, and indices on outgoing lines are paired with momenta as p^ ^ q'^ r^; 
Ai^ originates from the term where g^^ connects the e~-i/e, e^-i>e fermion lines. 



The soft part is now clearly separated from the remaining non-IR part, used in the CEEX 
exponentiation for construction of 0{a) corrections. We have ordered the expression, 
with the help of expansion similar to the contact interaction for W propagator as well. In 
05 [^] we use an auxiliary fixed transfer to, independent of the place where the photon is 
attached to the fermion line. In fact to is arbitrary and the choice to = could be used as 
well^. With the help of ^ we provide the residual contribution calculated as a difference 
of the expression calculated with the true t-transfers {ta or tb) and the auxiliary to one. 
Note that ^ = ^ + ^. Each of the contribution to the sum given in the first equation of 
(jHl) is independently gauge invariant. 

We can see, that it was possible to separate the complete spin amplitude for the process 
e~^e~ — »• z/eZ/e7 into six individually QED gauge invariant parts. This conclusion is rather 
straightforward to check, replacing photon polarization vector with its four-momentum. 
Each of the obtained parts has rather well defined physical interpretation. It is also easy 
to verify that the gauge invariance of each part can be easily preserved to the case of the 
extrapolation, when because of additional photons, condition Pa+Pb = Pc + Pd + ^i is not 
valid. Let us elaborate on this point a bit more. 

• The first term of the soft photon type So-/(fci)^ [^] is gauge invariant thanks to invariance 
of the standard ISR soft factor 3^-1(^1). It is also of the universal form, identical for the 
diagrams with s-channel Z-exchange as well as t-channel W. 

• The next two terms {rfji + A4^) and {rfn + M.^ ) originate only from diagrams of 
t-channel W exchange. For the gauge invariance to hold, the t-channel transfers have 
to be ta = {pa - Pc)^ h = {Pa - h - Pc^ for the first term (and ta = {pb - h- pdf, 
tb = {Pb — PdY for the second one)^. 

• The consecutive two terms rd, (^^0 ^^^ ^fn (actO ^^^ again of the same universal 
form as for any s-channel process and gauge invariant by construction. These are also 
the terms which lead to leading-log (but not infrared) singular terms, after phase space 
integration. In the corresponding Feynman diagrams the photon polarization vector and 
its momentum stand side by side. 

• Finally, for the last expression (7W^ + TW'^) to be gauge invariant it is enough that for 
the two terms choices for transfers ta, tb are identical; the same reduction procedure^ is 
used. 

Simplest case of e~^e~ -^ VfiV^ 

Let us finish this Section with the discussion of the Z exchange part of the amplitude 
(121) in simple language of spinors and four- vectors. This part of the amplitude is important 

^The choice is nonetheless important from the point of view of efhciency, it affects size of corrections 
in CEEX expansion scheme. Condition that limfcj_>o ta/b = ^0 is desirable. 

^It is interesting to realize that only part of the diagram is involved in the cancellation, namely fermion 
and boson lines, from which emission of the photon takes place. This observation will become useful in 
case of double bremsstrahlung amplitudes. 

^Mechanism of gauge cancellation is fulfilled already at the level of this part of bosonic interaction 
alone. This observation will be useful in study of double bremsstrahlung amplitudes. 



because it will define the framework for our main results collected in Section 3. 

"7 U (9) 

+eQe v{pb,Xb) C^ikl) ^^ M|/} u{pa,Xa) 

the superscript ac or bd can be dropped in M|/} as it does not depend on the t transfer 
of W propagator, then M{/} = TZz of (j211), see later in the text. 

The gauge invariance of the two sub-parts proportional to /(fi is straightforward to see, 
because the expression ^i /*^{ki) alone is gauge invariant thanks to f^if^i = 0. These 
parts of the amplitude do not contribute to infrared singularity, however do contribute 
to the big logarithm related to coUinear singularity (once amplitudes are squared and 
integrated over the phase space). That is why, we will call these parts of the amplitude 
as infrared finite coUinear singular. The remaining part of the amplitude: 
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K~i} (a'O =^Q^ ^(p^' ^^) ^\n r^ 4(^i) "(P- ^«) 

+eQe viPb, h) 4(A;i) ^^^ Mfjy n(p., A„) 
factorizes thanks to the orthogonality for Dirac spinors: 
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into gauge invariant soft photon factor and Born amplitude: 

''-b^^{ki,pa)Sx\a + T-r^KAki,pb)5x\^, 



2kiPa " 2kiPb 

(12) 



The gauge invariance, takes place in case of Z exchange (and also W exchange if 
approximation of contact interaction is used) because then Mj/} = ^Yn ~ ^m^ ^^^° 

two parts of sj-i (A^i) are diagonal respectively in indices AAa^b. The Born level spin am- 
plitude factorizes out, and the gauge dependent soft factors for emission from electron 
and positron lines, can be summed to gauge invariant si-j (/ci). For the explicit definition 
of bcr-^{ki,pb), see e.g. formula (231) of ref. [8]. This part of the amplitude is infrared 
singular. We will use the factorization of the soft factors explained here, also later in the 
paper. 



Note that in case of single photon Z exchange amphtude, we got only three gauge 
invariant parts: infrared-singular, and two others contributing to coUinear-singular terms 
(after phase space integration). The residual (contributing only non enhanced terms after 
phase space integration) terms are absent. Note however, that they are present in case of 
the W exchange diagrams. 

Finally let us comment that similar pattern of amplitude separation into gauge invari- 
ant parts can be observed for W^ -^ Ivf) in [13]. 

3 Double bremsstrahlung 

In the present section we will study the amplitudes for double bremsstrahlung in e^e~ — > 
v^Ve production process. There are two classes of diagrams in this case, the first one with 
Z boson exchange in s channel and the second one with t channel W exchange. Similarly 
as in previous section and single bremsstrahlung we will look if gauge invariant parts of 
the complete amplitude can be defined. We will be also interested, if this can be done in 
a semi-automatic way, directly from the Feynman rules. 

The presentation of the complete amplitudes is rather difficult because of their length. 
To avoid lengthy formulas, we will start with largely incomplete set of diagrams, nonethe- 
less sufficient to localize some gauge invariant group of terms. Once localized, it will be 
hidden under symbol l}^ and, to the remaining part contributions from the next diagrams 
will be added. Again gauge invariant group of terms will be searched for. This procedure 
will be repeated until the complete list of diagrams of our process will be exhausted. The 
choice for the first diagram in this procedure is motivated by its particular (unique) form, 
later it is motivated by the form of the gauge dependent rest remaining from the previous 
step. For short hand notations we will use extended subscripts and superscripts for L\. 
For example we will use symbol L \l ^(n), to denote contribution for the diagram with: 
first photon of momentum fci and second k^^ attached to incoming e~ line. The number 
n in bracket (if present) will denote that it is only a part of the contribution from the 
particular Feynman diagram (or diagrams). There will be often bar over this number to 
point that the particular part is gauge dependent. 

Let us start our iteration with diagrams involving double fermion propagator, that 
is diagrams where two photons are attached either to incoming electron or to incoming 
positron. These are the only diagrams with k\ ■ k^ term in fermion propagators. Our 
first aim will be thus to localize the parts which are gauge invariant by themselves, and 
include this k\ ■ k^ term. Let us consider the eight diagrams with the photon lines attached 
either to electron or positron line, see fig El Explicitly, we will write down the part of the 
amplitude corresponding to the incoming electron line only. In fact the diagrams with Z 
and W exchange are quite similar: 





Figure 2: Double emission from electron 
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Expression TZb describing final state neutrino interaction and Z or W exchange, will 
be defined later, respectively in formulas (j24j) and (j27|) . We can separate formula (J13p 
into the following parts 



^fci.fe _ i^iM^^-^ ^ L^i''=^(2) + lI]_'^'{3) + L^i'^^(4) 
where: 



(14) 



2k,p.-2k2P.-2k,k2 ^-^''^ ^ ^-('^^ "^^-'"^ 



+ -2k,p.-2k2P.-2k.k2 ^-^'^^ ^^^ ^-^'^) "^^-'" 

(15) 

is gauge invariant by construction, thanks to terms f^i ^„^{ki) and }i2 ^2(^2) which are 
gauge invariant by themselves. This is similar to the case of single bremsstrahlung. The 



second part 
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(16) 

is also gauge invariant, but it need to be checked by direct calculation. This contribution, 
like the previous one, is free of infrared singularity. In definition of LJ_' ^(2) we had to 
introduce subtraction; terms proportional to }i2{—^^ — ) and (j5a + rn){^—^:^ — ^^g^ — ) 
The subtraction terms are added back to (jT^ as formulas (fTTj) and (fTHj) . but with the 
opposite sign of course. It is important to realize, that these subtraction terms are defined 
uniquely by the Z exchange part of the amplitude for single photon emission^, see formula 
(P). It has to be multiplied by the soft photon factor for the second photon and separated 
into ifrared finite and infrared singular parts, as explained in section 2. 
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The next term L^_l ^(4) is, also free of k\k2 and its numerator is linear in the photon 
momentum 



^e-'^(4) = {eQefv{p,,\,)nB[-^^ 4(fc2) ^^ 4(A;i) U{j>a.\a) 
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^The subtraction term for the W exchange differs only by the replacement TZb = TZw- 
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The complete contribution from diagrams of fig. El formula lfT^ . is not gauge invariant. 
The last two terms are gauge dependent, but are also relatively short. The first one (fTTjl 
has a structure of Born amplitude multiplied by soft photon factors, the second one (fTH|) . 
has a structure of soft photon emission for one of the two photons only, see discussion at 
the end of Section 2. 

Once we have completed the diagrams with two photon lines attached to one fermion 
line, let us turn to another group of diagrams, where one of the photons is attached to 
electron and another one to positron line, see Fig. IHl Note that for subgroup of diagrams 
with Z boson exchange only, these are the last contributing diagrams: 
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(19) 





Figure 3: Single emission from electron and positron 
As in the previous case the expression for i^„i' + can be easily separated into parts. 



^k,_M ^ L^i-^^^41) + l';j%{2) + Lll^'%{3) (20) 
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The first part LJ_' ^+{1) is by construction gauge invariant, it is also the only part from 
this group of diagrams with numerator proportional both to the momenta of ki and ^2^ 
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The second term LJ_' ^+{2) has two contributions, of numerators hnear either in ki or A;2, 
it reads: 



-2kipa 



-2k2Pb 


^^ 7^ 


-2A;ij9b 


^^ 71 
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- ]^b + m 



(23) 



(22) 

+v{Pb.\) fa.ih)^^ ^^^^S ^-^^^) ''^^"'^'^) 

+v{Pb,Xb) c(^i)^^ ^^^i£: ^^^^'^ ''^^"'^'^^ 

Finally the third one LJ_' ^+ (3) is free of the fci or ^2 in the numerator. 

ie-,t+(3) = (eQe)^(t'(p., A,) ak2)^^^ ^^^^ ^-^^^^ "^^'^'^'^^ 

+.XP..A.) 4(*i)^^;;^ '^B^ 4(« "(Pa, A.) 

To complete the sub-set of diagrams for double bremsstrahlung from initial state, the 
contribution of the double emission from positron line should be added. We will omit ex- 
phcit formulas here, the explicit expressions for L\^ ^(1), L\^ ^(2), Lg+ ^(3), L^^ ^(4) can 
be obtained from L^i' ^(1), L^i' ^(2, )i^gi' ^(3), L^J ^(4) by analogy or explicit calculation. 

3.1 Diagrams with Z exchange 

Before going to the more complex case of W exchange, where complications due to t 
dependence of W propagator occur, let us concentrate on the simpler one, Z exchange 
alone. The diagrams discussed so far, represent then the complete gauge invariant am- 
plitude for the process e^e~ -^ 1/^1/^77. In such a sub-group of diagrams (for the process 
e^e~ — > Vf,Vf,ll) symbol 71b=z represents always 

nz= [rivl + a^''))^ {u{p,,K)^^{vl + a^'')v{pdAd))BWz{{Pc + Pdf) (24) 

which is a constant algebraic expression, independent on photon momenta and identical 
for all diagrams listed. The Z boson propagator BWz^ipc+VdY) depends on the invariant 
mass of the outcoming neutrinos only. The bi-spinor indices of 7^^, 7^7^ matrices which 
enter into the matrix products of formulae such as (fTT^ to (^Hjl are explicitly given as a/3. 
The complete amplitude reads: 

e~ e+ e~,e+ 

=Lll!^\l) + L^i''^2(2) + L^i''=^(3) + Ll'/\l) + 
L'l^/^l) + i:^v''(2) + i:'l''^(3) + L^V'^(4) + 
L^^J%{l) + L';j%i2)+L';j%{?,). 

12 



Where L i' ^ is given by formula (jl3|) (or if separated into parts by (J14j) ) and L i' + by 

(fT^ (or (j2ni))- For L^l' ^ the expressions of L^J ^ can be used, with appropriate changes 
of signs momenta etc. 

The formula for the complete spin amplitude {Z exchange only) can be easily re- 
ordered into consecutive contributions A^i, A^2,-A^3, •••, each gauge invariant, and each 
expressed as group of L's in the bracket or just individual L: 

=Mi + M2 + M3 + Mi + M5 + M6 + My 

=L^i''^^(l) + ill'"' (2) + LJV''(1) + ^ei'''(2) + ^e-t+(l) ^^^^ 

+ {Lll''\A) + Lll''\A) + L^i't+(2)) + ((^e-''(3) + ^'^''(3) + ^'-'+(3)) 

As one can see, the sum of terms Mi to M5, contributing to /3| of the CEEX exponen- 
tiation scheme (these terms are not infrared singular at all), is gauge invariant and clearly 
separated from the rest. It can be even further sliced into five parts, each individually 
gauge invariant. The last two terms, Mq and My corresponds respectively to Pi and Pq 
(multiplied by one or two soft photon factors) and can be obtained from lower order of 
perturbation expansion. It is rather straigtforward to see, that the term My consist of 
Born level amplitude multiplied by soft factors corresponding to emission of two photons. 
The term Me can be seen as consisting of two factors; for one of the photons soft factor, 
and for the other one term of Pi- To see it better it is convenient to order the expression 
accordingly to terms proportional either to f^i or ^2- 

Note also, that for the each of the parts to be gauge invariant, it is not necessary that 
four-momentum conservation is fulfiled. That is why, the separation is easily adaptable 
to extrapolation procedure as used in KKMC. 

This completes our discussion of results for s-channel exchange of Z. Let us now turn 
to the contributions related to the t-channel P^-exchange. 

3.2 Diagrams with W exchange 

First, let us note that all formulae presented so far are valid for the diagrams involving W 
exchange as well. The difference is that instead of formula P^ for TZb one should use: 

TZw = (7m(1 - I'HPd, Ad))^ (uip,, Ae)7^(l - 7')) BWw{t) (27) 

The spinorial form of this expression is universal, and as in the case of Z exchange, in 
all places the same expression is to be used. The difference hes in t dependence of W 
propagator, the transfer will depend on the way how the photon lines are attached to the 
fermion ones. Nonetheless in some groups of terms gauge cancellation do occur anyway 
in the same way as before. If we recall the part of the W exchange amplitude, written in 
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analogy to pH|l as: 



M 



A 



--Ml + M2 + Ms + M4 + M5 + Me + M-, 



' ki ,k2 



-ki,k2 



L'^].'"'{1) + L'll:'\2) + L^r^(l) + L^l''^^(2) + L^i'7+(1) 



ki,k2i 



-k\,k2 



e ,e^ 



(28) 



'k\,k2 , 



'k\,k2 



k\,k2 I 



'k\,k2 I 



'ki,k2 



+ {L'^l''\A) + L'f\l) + l11:':^{2)) + {{L':j'\1) + ^^^''=^(3) + Ll]:%{l)) 



then the parts Mi, M2, M3, M4 and M5 remain gauge invariant. Only the last 
two terms will need contributions from diagrams with triple and quatric gauge boson 
couplings for the gauge invariance to hold. To visualize this point the bar sign is now 
placed over Mq and M7 . Note, that as already pointed in the previous subsection these 
are the contributions, which could be obtained (extra soft photon factors are only needed) 
from the results of the calculation at lower perturbative order if the complications due to 
variation of W exchange transfers were not taken into account. 





Figure 4: Single and double emission from W 





Figure 5: Four boson coupling and coupling for unphysical x field. 
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Let us continue with the second part of our discusion now. The two (now gauge 
dependent), contributions Mq and A^y will be first completed with diagrams from the 
left-hand side of fig. HJ Note that these diagrams are the last ones with photon line 
attached to incoming electron/positron, thus the last ones contributing collinear and/or 
soft singularities. The contribution to the scattering amplitude from these new diagrams 
reads: 

L'e-'w ={eQe?BWw{{Pc + k2- Paf)BWw{iPc + h + h- paf) 
v{pb, Xb) (1 - 7^)7'' v{Pd, Xd) 
[9^.AP - (l)p + 9up{.q - h)^ + g^pih - p)u] (4i (^1))'' (29) 

uip,, Ae) (1 - 7^)7^^ ^^l^~f' 4(^2) uip,, Aj) 
+ (1^2) 

Here p = p^ - p^, = -{p^ _ p^ + fc^ + k2) and q = p^ - p^ + k2 = -{pd - Pb + h)- 
Similarily one can write contribution LJ^^ for the two diagrams with emission from 

positron and W ^ we will omit the corresponding formulae. As before, we separate L^i' ^ = 
LJ_'^(k^) + LJ_'^{k^) into parts, (k^) marks contribution where only f^2 is taken from 
the fermionic propagator and (A;") marks the rest. The explicit formulas are: 

Ll-'wik^) HeQeYBWwiiPc + k2 - PaY)BWw{{Pc + ^2 + A;i - p^Y) 
v{pb,Xb) {1-1^)Y v{pd,Xd) 

[g^uip - q)p + 9up{q - h)^ + g^pik^ - p)^] {el^{ki)Y (30) 

«(p„ A,) (1 - 7^)7^^ ^^ f.,ik2) uipa, Aj) 
+ (1^2) 
and 

L'-'wik') HeQefBWwiiPc + k2 - PaY)BWw{iPc + k2 + k, - Paf) 
'v{pb,Xb) (1-7^)7'' v{pd,Xd) 
[9pu{p - q)p + gvp{q - ki)^ + g^piki - p)^] {el^{ki)Y (31) 

u{p,, A,) (1 - 7^)7^^ -=^ 4(^2) «(Pa, A,; 
+ (1^2) 
Let us start with the second one, which can be easily transformed (with help of Dirac 
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equation) into 

Lt-%{k') ={eQe)'BWw{{Pc + h- Paf)BWw{{Pc + h + h- p^)') 
(v{pb,Xb) {l-l^)Y v{pd,Xd) 
[g^iuiPd -Pb-Pc + Pa- k2)p + 9upi-'2h)f, + gf,p{2ki - p«)^] (e^^ {ki)y (32) 

u{p,, K) (1 - j')Y 3^ 4(^2) u{p,, A,)) 

+ (1^2). 
This contribution can be separated even further: 

where 

Lt-'wi^) ={eQeYBWw{{Pc + k2-Paf)BWw{{Pc + k2 + h -Paf) 

{v{pb, A,) (1 - 7^)7'^ v{p,, A,) u{p,, Ae) (1 - 7^)7'^ 3=^ 4 (^2) «(Pa, A«)) 

(Pd -Pb-Pc + Pa- k2) ■ e*^ (/Ci) 
+ (l--2), 

(34) 

and the second term 

^e-'l^(2) =(ege)'5iyH'((Pc + k2- Paf)BWw{{Pc + k2 + h - Paf) 

y{pb,Xb) (1-7^)7'' v{pd,Xd) 

V-(C(^i)).(^i)m+(C(^i)),(^i)-] (35) 

«(j9„ A,) (1 - j')Y ^^ CM u{p,, Aj) 



-'2k2Pa 

;i-2) 



is gauge invariant by itself. The third one is exphcitelly less divergent in collinear config- 
uration: 



^e-t^(3) = - {eQ,fBWw{{Pc + k2- Pa?)BWw{{Pc + k2 + k, - Pa)') 
'v{Pb,Xh) (1-7') C{ki) v{pd,Xd) 



M(Pe, Ae) (1 - 7')^a ^7^— 4(^:2) «(Pa, A„) 
—ZK2Pa 

+ (1^2). 
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Let us now turn to the other part of the amphtude, and again present it in a form of the 
sum: 

^e-'^(^°) = L'l:%{l) + Lji;^^(5) + Lji;^^(6) (37) 

where first term reads 

- fcl,fc2 



'^a + 1 
~'2k2Pa 



v{p,, A,) (1 - 7^)7^ v{p,, A,) u{p,, A.) (1 - 7^)7^ ^^ 4(^2) «(Pa, Aj) 



{Pd -Pb-Pc + Pa- h) ■ 4i (^1) 

+ (l--2), 



(3^ 



the second one 

Ll'j%{h) ={eQefBWw{{Pc + k2- Paf)BWw{{Pc + k2 + h- pa)^) 
(v{pi>,Xb) {1-i^h^ v{pd,Xd) 
2[-(4(^i)).((^i)m - iPb),) + {el.ikr)) ^{k^l] (39) 

u{p,, A.) (1 - 7^)7'^ ^g-^ 4(A;2) u{p,, Aj) 
+ (l--2), 
and the third 

Le-'wi^) ={eQeYBWw{{Pc + k2- Pay)BWw{{Pc + k2 + k,- PaY) 
'v{pb,Xb) (1-7') 4(A;i) v{pd,Xd) 

-2k2Pa 

+ (1^2). 



u{p,, A,) (1 - 7')(^2 - ^a) ^^ 4(^2) u{p,, A„)) 



(40) 



The last two terms can be modified further, and some terms neglected. One can check 
that these terms contribute at the level of ^ only, and that is why, we will exclude them 
from explicit considerations. After these simplifications, we finally obtain gauge invariant 
by itself 

L';'-%{5) ={eQefBWw{{Pc + k2 - Pay)BWw{iPc + k2 + k,- p^Y) 

[v{pbAb) (1-7^)7'' v{pdAd) 

2[-«(^i)).(^i)m+(4(^i)),(^i).] (41) 

«(Pe, A,) (1 - 7^)7^^ ^g-^ 4(A;2) «(p.. A,)) 

+ (1^2). 
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and now explicitelly less divergent in coUinear configuration term 

^e-'^(6) ={eQefBWw{{Pc + k2- Pay)BWw{{Pc + k2 + h- Paf) 

'v{pi„Xb) (1-7') Cih) v{pdAd) 



U{p,, Ae) (1 - l')h ^g^ f.^ih) Uipa, Xa)) 

+ (1^2). 

Let us now turn to the diagrams of double emission from tlie W, see fig. 01 riglit-liand 
side. Tlie corresponding amplitude can be written as: 

L^;^ ={eQefBWw{{Pc-Pa)')BWw{{Pc + h - p^f) 

BWw{{Pc + ki + k2-paf) 
(v{pi>, Ab) (1 - 7^)7^^ V{pd, \d) u{Pc, Ac) (1 - 7^)7" U{pa, Xai) 

[gauij) - q)p + Qupiq - ki)^ + g^rpiki - p)^] (^3) 

[<(p' - q')p' + gl'iq' - ^^2)^ + g^p'ik^ - p'Y] 

+ (1^2) 

where p = -q' = Pd - Pb + k2 = -{pc - Pa + ki), q = Pc - Pa = -{Pd - Pb + ki + ^2) and 
P' = Pd - Pb = -{Pc -Pa + ki + k2). 

As usual, we will represent this expression in a form of the sum: 

^&:S = ^&:5?(i) + ^&:5?(2) + ^S-.-^is) + ls.:S(4) + ls.;S^(5) + lJ^;^(6). (44) 

The first term, proportional to Born amplitude multiplied by factor depending on polar- 
ization of the two photons takes the form: 

^feS(l) HeQefBWw{{Pc-Paf)BWw{{Pc + k^-paf) 
BWw{{Pc + ki + k2-paf) 

%{pb, \b) (1 - 7')7'' v{pd, \d) u{pc, \c) (1 - 7')7/. u{Pa, Xa)) (45) 

\p-q)-<Aki) {p' - q') ■ {CAk2)) 

+ (1^2). 
Terms where dependence on polarization of only one photon factorizes out from the 



amplitude take a form 



L';^Pi2)=ieQ,)'BWw{iPc-Paf)BWw{iPc + h-pay) 
BWw{{Pc + ki + k2-Paf) 

(v{pb, Xb) (1 - l'')Y v{Pd, Ad) Uip„ A,) (1 - 7^)7" U{pa, Xa)) ,^g^ 

ip' - q') ■ (et^ih)) 
and 

Lw^i^) HeQe^EWwiiPc- Pay)BWw{iPc + h-paf) 

BWw{{Pc + ki + k2-Paf) 

{v{pk, Xb) (1 - 7^)7'' v{pd, Xd) u{p^, AJ (1 - j^)Y u{pa, A j) 
2[-{e:,{k2)Uk2), +ik2UCAk2)),] 

ip-q)-<A^i) 

+ (1^2). 

The two last terms ()46|47p are partially gauge independent, respectively for the polariza- 
tion vector of the first and second photon. The remaining fully gauge dependent parts of 
the amplitude read: 



L'^fiiA) = - ieQ,yBWw{iPc-Pay)BWw{iPc + k.-p^f) 
BWw{{Pc + ki + k2-Paf) 
(v{pb, Xb) (1 - 7^) h v{pd, Xd) u{p,, X,) (1 - 7^) ^ 



^Ww{{Pc + ki + k2-Pa)'} 

v{pb, Xb) (1 - 7') h v{pd, Xd) u{p,, AJ (1 - 7') 4 (k,) u{pa, A J ] (4S 
P' - q') ■ {C{k2)) 



' ip' - q') ■ iCAk2)) 
1^2), 



Lw'wi^) =ieQe)'BWw{iPc-Pa)')BWw{iPc + ki-pa)') 

BWw{{Pc + ki + k2~Pa)^) 

{viPb, Xb) (1-7') 4 (^2) vipd, Xd) u{p,, A,) (1 - 7') ^1 u{pa, Xa)) (49) 

{p ~ q) ■ ^A'^i) 
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Finally 



L'^wiQ)=ieQefBWw{{Pc-Pa)')BWw{{Pc + h-pay) 
BWw{{Pc + ki + k2-pay) 

(v{pb, Xh) (1 - I'h'" v{Pd, Xd) U{pc, Xc) (1 - 7')7' u{Pa, Aa)) 

[gup{q-ki)^ + g^p{ki-p)^] (e^,(A;i))^ 



+ (1^2) 



At the last step, let us turn to contributions from diagrams presented in fig. 13 The 
diagram with contribution from quatric gauge coupling reads: 

L%'t ={(^QefBWw{{Pc- Pa?)BWw{{Pc + k2 + ki -paf) 

(^(pb,A6) (1-7^) 4(A;i) v{pd,Xd) m(p„ Ae) (1 - 7') f^,{k2) u{pa, Xa) 
+v(p6,Ab) (1-7^) 4(A;2) v{pd,Xd) u{pc, Xc) {1 - 1^) f^^{ki) u{paAa) 

+v{pk, Ab) (1 - 7^)7^^ v{pd, Xd) u{pc, Ac) (1 - 7') 7m u{pa, Xa)el_^ (fci) ■ e^^ (/cs)) . 

(51) 

It is convenient to write it as a sum of two parts 

^^'^ = ^^'^(i) + ^'^''^(2) (52) 

where 

L'^^(l) =2{eQefBWw{{Pc-Paf)BWw{{Pc + k2 + k,-paf) e^^^k,) ■ el^{k2) ^^^^ 

V{pb, Xb) (1 - 7^)7^ V{pd, Xd) U{pc, Ac) (1 - 7^) 7m ^(Pa, Xa) 

and 

L^','^(2) ={eQe)^BWw{{Pc~Pay)BWw{{Pc + ^2 + ki - Pa)^) 

(^V{pb, Xb) {1 - 1') CC^l) v{pd,Xd) U(pe, Ae) (1 - 7') f^,{k2)u{pa,Xa) 

+^(pfc,Afe) (1-7^) 4(A;2) v{pd,Xd) u{pc, Xc) {1 - 1^) f„^{ki) u{pa, Xa))- 

( 

Contribution from the diagram involving internal x li^is reads: 

L'^i' =(eQe)'M^ BWw{iPc-Pa)')BWw{{Pc + ^2 + k, - p^)') 
(v{Pb,Xb) (1-7') C(^i) ^iPd,^d) 

U{p„ Ae) (1 - 7') ^2(^2) U{pa, Aj) 

+ (1^2). 



(54) 



(55) 
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This closes the hst of all diagrams entering the complete spin amplitude for the process 
e^e~ — >■ z/eZ/e77. The contributing terms were obtained from the Feynman rules, and were 
grouped on the basis of rather straightforward rules; gauge symmetry and nature of 
singularities in infrared and coUinear limits (phase space integration was not necessary). 

The complete gauge invariant part of the spin amplitude of W exchange can be now 
written as: 

My, = M^ + Ml,, (56) 

where M.y^ (technically identical to the amplitude of Z exchange) was given by formula 
(1^ . and new part, specific to the W bosonic interactions, reads: 

Tkl,k2/T\ I rfcl,fc2/oN I 7-^1,^2 

We can now write the complete spin amplitude, of W interactions, as a sum of gauge 
invariant parts: 

Ms + Mg + Mio + Mil, 

where 

Ml =L';j''{l) 
M2 =L^'i''^(2) 
Ms=L';/\l) 
M,=L';/\2) 

M, =L'^^'\l) + L'f\l) + Ll^_:%{2) + Ll^%{l) + LlX%{l) 

Mr =L'^J'\^) + L^i''=^(3) + Ll^J%{^) + lI]:%{1) + L^\^^(4) + L%''^{1) + L^'^(I) 

Ms=L\ 



k\,k2 



M, =L':X%i2) 
Mio =L':i%{5) + L^^ib) + Lf^;^(2) + L^_$(3) 
Mil =L':l.%{3) + L^X^iS) + L^^i'5.(6) + L^m + L^^;^(4) + Lg^;^(5) + L^;J^(6)+ 

T ki .ko /ri\ . T k-i .ko 



Tki,k2/r)\ I Tki,k2 



'w,x ■ 

(59) 



21 



The matching of the L'^{n) terms into gauge invariant parts Aii of the amphtude 
is straigtforward and based on the type of singularities present/absent in the particular 
group. Each of the listed below contributions TWi-TWn can be given some physical 
interpretation. In some cases, appearance of such parts may seem rather unexpected. 
In brackets we provide symbols such as (lA), they denote the name of variables used in 
KKMC [7] Monte Carlo, as keys for the parts of the amplitude: 

• TWi (lA) , contribution of the infrared non-singular contributions of double emission 
from electron line, part with straightforward gauge cancellation within the terms 
originating from diagram of two photons attached to the same incoming electron 
line. 

• A^2 (IV2) contribution of the infrared non-singular contributions of double emis- 
sion from electron line, part with non-straightforward gauge cancellation within the 
terms originating from diagram of two photons attached to the same incoming elec- 
tron line. Part of the diagram contribution had to be subtracted; more precisely 
expression without kik2 product in electron propagator. This subtraction term is 
recupered in A^g and Aij. 

• Ai^ (Ik), Aii (IVl) as in previous two cases but for emission from positron line. 

• 7VI5 (18) infrared non-singular contributions of single emission from electron- and 
another single emission from positron line. This contribution is gauge invariant by 
construction. 

• A4q (I9X) , (I9Y) , (I9Z) , (I9T), part of the amplitude with infrared factor for 
one photon, and for the second one infrared non-singular gauge invariant contribu- 
tion. For the diagrams with W exchange contribution from diagram with photon 
emission from W need to be taken. For the gauge cancellation to hold, relation 
between t-channel transfers in W propagators and momenta multiplying photon 
polarization vector need to be fulfiled. Nonetheless certain freedom in choice is left. 
It was useful in construction of extrapolation procedures^. 

• Aij (IVI) part of the amplitude with infrared factors for both photons. For the 
diagrams with W exchange contribution from diagrams with single and double emis- 
sion of photons from W needs to be taken, also part of the diagram with quartic 
gauge coupling was needed here. 

• Ais (171) , (172) part of the amplitude with infrared non-singular contribution of 
emission from electron for one photon and for another one part of emission from W 
which is self gauge-conserving. 

• Aig (171) , (172) as in previous case but for emission from positron. 



*^Idcntical condition, also originating directly from Ward- identities, need to be preserved in Alio and 
the similar one in Mn- 
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• A^io (I9sl) , (I9s2) part of the amplitude with infrared factor for one photon 
and for another one part of emission from W which is self gauge-conserving part. 

• A^ii (19) , (I9B) , (110) all remaining parts, they turn out to be free of singular- 
ities both in collinear and soft limits. 

• Let us comment that in the limit Mw -^ c>o all contributions from Ai^ to Ain 
disappear. In this limit amplitudes for s-channel Z exchange and t-channel W 
nearly coincide. The only remaining difference is the coupling constants and hard 
interaction part of the amplitude given respectively by formulas ()24j) and (P7j) . This 
is an extension of similar observation of reference [14] instrumental in construction 
of extrapolation procedures of ref . [9] , to the case beyond real photon interactions 
with fermions only. 

• Let us point that in many places we have used separation of the WlVj vertex into 
three parts; (i) the one with the g^u tensor along line connecting fermion lines, (ii) 
the part internally preserving gauge symmetry, (iii) the remaining part which we 
often could reduce significantly with the help of Dirac equation (for the fermion 
lines connected with the WW'j vertex by W propagator. 

• Finally let us note, that the above separation into gauge invariant parts can be 
continued even further. For example it is rather easy to separate A4q into four 
parts. For each, emissions of individual photons are attributed either to electron or 
positron line. 

Let us note that we have not exploited to the end the properties of A^u. It was not 
interesting from the point of view of our main purpose, which is implementation of the 
matrix element to the environment of Coherent Exlusive Exponentiation. Also in case of 
At 11, contrary to the cases A^i to A^io, similarities with first order results could not be 
seen. This is rather natural, as for example quatric gauge couplings are absent in first 
order. In this case hint on pattern of constructing amplitudes of even higher order using 
iteration techniques could not be found. To this end, discussion of the amphtudes of 
triple photon emission would be needed. If conclusive, it would point to solutions beyond 
next-to-leading-log approximation, thus beyond imminent interest of the present paper. 

The gauge invariance was not the only element of the criterium which was used here 
to split amplitude into gauge invariant parts. Equally important was that the two main 
sources of the radiation, incoming beams, form the unambiguous frame with respect to 
which, photon energy and the angles of photons with respect to fermions could be defined. 
That is why there was no need to make any reference to the regulators, singular terms 
could be localized already at the amplitude level and in fully differential manner, with 
no need to partially integrate phase space. The expansion in the contact interaction for 
W propagator enabled to place the gauge cancellation effects of emission from t-channel 
W within the frame of ISR radiation. Also relation between amplitude for double and 
single photon emission had to be exploited to close down window for ambiguities. Once 
these assumptions and properties were exploited, the solution seem to be unique, up to 
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may be grouping or further splitting of the obtained parts. Confirmation, whether this 
is accidental property and observation which hold for the particular cases and up to the 
second order only, may require calculation to be extended to at least third order. 

4 Some points on extrapolation 

Let us summarize here some specific issues related to extrapolation procedure of CEEX 
scheme described in detail in ref. [8] for purely s-channel hard process. One of the im- 
portant property of perturbation expansion, rearranged to improve convergence into ex- 
clusive exponentiation, is that parts of the amplitudes need to be appropriately shifted 
between the rearranged orders of expansion. We will concentrate on issues related to real 
bremsstrahlung only. In particular parts of the higher order terms (directly calculated in a 
standard way) which are already available at lower level of CEEX perturbation expansion 
need to be localized and subtracted in a clear way. Only remaining residual parts, called 
f3^, f3^, (3"^ etc. [2] will be indeed the term of the given newly rearranged order. The use 
of 13 functions is unambiguous if sufficiently high order of standard perturbative calcula- 
tion is available to calculate matrix element, for the configuration with all real photons. 
However it is not always the case, practical solutions for exponentiation require definition 
of methods how to calculate matrix elements for the kinematical configuration with large 
number of real photons, using results of first (or second) order of perturbation expansion 
only. 

There are several rules which extrapolation procedure must fulfill. Already the lowest 
order must include all terms with the highest power of infrared singularity and for all 
kinematical configurations of arbitrary number of real photons and in a fully exclusive 
manner. Then, first order provide all terms with next to highest power of infrared singu- 
larity, etc. Let us stress that reduction/extrapolation procedure of exponentiation offers 
some freedom of choice. This freedom can be used to further improve convergence of 
perturbation expansion. The best guidance is of course comparison with result of even 
higher order of expantion to minimize their contribution. If such results are not available, 
higher order leading log results can be used instead. Finally, let us stress that if suffi- 
ciently high order of perturbation expansion is available, dependence on particular choice 
of extrapolation drops out and unique result, identical to the one of direct perturbation 
expansion without any reordering, will be obtained. Unfortunately this is not expected 
to be the case in foreseeable future. 

In case of diagrams with Z exchange the question of choice of extrapolation procedure 
is straightforward. Inspection of first order (formula EJ and second order (formula I^U]) 
amplitudes points to the following solution: the terms M.i to M.^ of (J2S1) should con- 
tribute to /3^, whereas the last two terms M.q and M.j can be directly obtained from the 
lower order. The M.^ can be obtained from (3^ by multiplication with the soft photon 
factor for the other photon. The (3^ can be identified as this part of M.^ (see formula 12)) 
which is proportional to /^i. The M.j can be obtained from the lowest order Born spin 
amplitude (3^ of Z exchange by multiplication with two soft photon factors, for each of 
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the bremsstrahlung photons, exactly as it should be in exponentiation prescription. The 
factorization properties can be easily seen if rather trivial manipulation on Dirac algebra 
is performed. 

In case of diagrams with W exchange, the question of choice of extrapolation procedure 
is slightly more complex, because of dependence of the transfers used in W propagators 
on photon momenta. That is also the reason why triple and quatric gauge couplings need 
to be included in the considerations. If the kinematical configurations of more than two 
explicit hard photons are taken, then the transfers calculated for W propagators can be 
defined in several ways. Our choice used at present in KKMC [7] is inspired by leading 
log considerations. For lowest order (/9°,) and if there was no addition photons, transfer 
to, can be calculated either as (i) to = {Pc — PaY or (ii) to = {Pd — PbY- If there is a photon 
coUinear to pb the first choice is closer to the transfer dominating higher order (i.e. single 
bremsstrahlung) spin amplitude. In general the choice (i) is thus more favored if total 
four-momentum carried out by the sum of all photons is pointing rather into direction of 
Pb than Pa- Otherwise the second choice is better. In case of single (or double) photon 
emission the choice how the transfers are calculated is basically the same. The only 
difference is, that the photons explicitely included in the particular contribution to /9^ 
or /5^, should contribute to the sum of photons mentioned above. The choice which pair 
of four momenta {pa, Pc or pb, Pd) is used in calculation for transfers, must be taken in 
calculation of algebraic expressions originating from direct W interaction with photons, 
for gauge invariance to hold. 

5 Summary 

We have presented complete results for the spin amplitudes of e'^e'^ -^ Uei^ell process. 
Using gauge transformation, as well as expansion with respect to contact approximation 
for W exchange, we were able to identify gauge invariant parts of amplitudes of the well 
defined physical properties. In particular the parts proportional to inverse of photon en- 
ergies (i.e. corresponding to infrared singularity), remaining parts proportional to inverse 
of the product of fermion and photon momenta (i.e. of the type of collinear singularity) 
as well as residual finite parts could be grouped together in a rather natural way. By 
comparison with amplitudes for the diagrams involving s-channel Z exchange we were 
able to observe certain pattern of universality for many of those terms. 

Let us stress, that some of the results presented here, could be expected from the 
properties of U(l) gauge symmetry and the corresponding Ward identities. They are 
known already since a long time, also in the context of QCD. The purpose of the present 
paper is mainly technical. We illustrate the scheme of step-by-step gauge cancellations 
and how they work for spin amplitude techniques. Finally, we show, how they helped to 
develop extrapolation procedures used in the KK Monte Carlo in case of neutrino channel. 

Let us point to the similar observation [15] as ours, in case of the single loop corrections, 
again for e+e^ -^ z/gZ/g process and also for e^e~ — > Ve^eH. 
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